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Abstract 

We study the fermion Dirac sea along the minimal energy path between 
topologically distinct vacua in the electroweak theory. We follow in detail the 
interplay between the bound state and the continuum in the fermion crea- 
tion/annihilation, as one passes from one vacuum to another. We calculate 
the quantum correction to the classical energy of the sphaleron barrier, arising 
from the fermion sea, and its contribution to the baryon number violation rate 
for non-zero temperatures. We find that the fermion fluctuations suppress that 
rate if the mass of the top quark is large enough. 
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1 Introduction 



The fact that the potential energy of the Yang-Mills fields is periodic in a cer- 
tain topological functional of the gauge fields, called the Chern-Simons number 
Ncs, was shown quite a long time ago by Faddeev and Jackiw and Rebbi 
In an unbroken Yang-Mills theory (like QCD) the classical energy barrier 
between the topologically distinct vacua can be made infinitely small due to the 
scale invariance. In a spontaneously broken theory (like the electroweak one) 
the scale invariance is explicitly violated by the Higgs v.e.v., and the height 
of the barrier is of the order of mw/a, where mw is the VF-boson mass and 
a = 5^/471 is the SU{2) gauge coupling. 

Transitions from one vacuum to a topologically distinct one over this barrier 
cause a change in the baryon and lepton number by one unit per fermion 
doublet due to the axial anomaly Q. Hence, this transition is a baryon and 
lepton number violating process and therefore of great physical significance. 
Although it is suppressed under ordinary conditions |^], its rate can become 
large at high densities [Q, |5|, high temperatures or high particle energies 
[0, ^. Thus, it could have occurred in the early universe [^, |l^, and an exact 
determination of the transition rate is important to understand the excess of 
matter over antimatter. 

The static configuration of the Yang-Mills (YM) and Higgs (H) field cor- 
responding to the top of the barrier and having Nqs = 1/2, called sphaleron, 



was first found by Dashen, Hasslacher and Neveu |11| and, in the context of the 
electroweak theory, rediscovered by Klinkhammer and Manton [12|; its energy 
depends somewhat on the H-to-W mass ratio. Moreover, a continuous set of 
static configurations of W- and //-fields, leading from a vacuum at Nqs = 
to another at Nqs = 1 and passing through the sphaleron at Nqs = 1/2, was 
presented by Akiba, Kikuchi and Yanagida (AKY) (see also ref. [|l4|) who 
minimized the Yang-Mills and Higgs potential energy for a given value of Nqs- 
This set of configurations describes a pass in the potential energy surface over 
the Hilbert space of fields, the sphaleron corresponding exactly to the saddle 
point of the pass. As most of the workers in the field, we shall simplify the 
actual electroweak theory by neglecting the Weinberg angle, i.e. reducing it 
to the pure SU{2) case (without the C/(l)) and considering the simplest case 
of one Higgs doublet. It should be mentioned that the sphaleron barrier has 



also been found for non-zero Weinberg angle |15]; however the difference to the 



2 



idealized case does not seem to be significant. 

The purpose of this paper is to study fermion field fluctuations around the 
whole sequence of static YM-H configurations leading from Ncs = to 1. 
To this end we investigate the behaviour of the Dirac sea of fermions as one 
passes the sphaleron barrier. This is of great methodological interest, since, 
as mentioned above, the baryon and lepton numbers are changed during this 
passage owing to the axial anomaly. We follow in detail this phenomenon. The 



discrete level has previously been studied in The complete Dirac spectrum, 
however, is calculated for the first time in this work. We apply these results to 
determine the fermionic contribution to the energy of the configurations on the 
path from Nqs = to Ncs = 1- This contribution is suppressed by the factor 
a = g^/47r ~ 0.04 relative to the classical one, but it also contains the large 
factor Nf which is 12 for the real world. It is reasonable to consider the range of 
parameters a <^ aNj <^ 1 where the fermionic energy is a first order correction 
to the classical energy, but might have a considerable numerical value. In the 
case aNf ~ 1 one would have to solve self-consistent equations for a sphaleron 
coupled to fermions. 

Having calculated the Dirac spectrum we evaluate the temperature depen- 
dent corrections of the baryon number violation rate due to the fermionic fluc- 
tuations. For a realistic set of parameters g = 0.67, A'^^ = 12 {2Nc -|- 3 = 9 
massless and Nc = 3 massive doublets), mn = mw, and mt between l.bmyy 
and 2.5 mw we find numerically that the fermionic contribution to the expo- 
nent of the Boltzmann factor increases the classical exponent by up to 65% in 
the temperature range O.Smw < T < 1.5mw, which can considerably sup- 
press the baryon number changing sphaleron transition rate. The suppression 
increases dramatically with the top quark mass; at mt ~ 3mw the fermionic 
contribution to the exponent becomes as large as the classical one. 

Finally our aim is to compute the baryon number and -density. We confirm 
the law B = Ncs and show that for Ncs ~^ 1 the created baryon becomes a 
free, delocalized particle. 

The paper is organized as follows: After recalling briefly the AKY [13| con- 
figurations along the barrier, we introduce the Dirac operator in the background 
field of those configurations and diagonalize it in a spherical basis. We show 
how the fermionic energy contributions can be extracted from the Dirac spec- 
trum. The energy of the Dirac sea is known to be divergent. Using a proper 
time regularization scheme, we separate the divergencies and combine them 
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with the classical energy whose parameters become renormalized. We next 
consider the high temperature transition rate and derive the fermionic part of 
its prefactor. The temperature range in which the result can be applied will 
be discussed. Finally we present our numerical results and discuss them in the 
last section of the paper. 

The simplified version of the standard electroweak theory that we are go- 
ing to study describes the SU (2) YM field interacting with the H doublet, 
and one left-handed fermion doublet showing a Yukawa coupling to two right- 
handed singlets. In principle, an SU{2) theory with an odd number of left- 
handed fermion doublets is not properly defined because of the so-called global 
SU{2) anomaly |17]. However, as we shall see in practical terms we do not 
encounter any problems related to that anomaly in the present work. 

The Lagrangian is thus: 

^ = C-YMH + i^F, 

J^YMH = -^F^,,F'^^^ + {D,<^>)\D^<^>)-^(^^^^-^^ , (1.1) 
Cf = i/jLh^D^i^L + XRij^d^XR-xPLMXR-XnM^L (1.2) 

with F^j^ = df^A"^ — dyA'^^ + £°'^'^A^^A1, and the covariant derivative being defined 
as = (9^ — iA'^^r"' /2. M is a 2 x 2 matrix composed of the Higgs field 
components $ = (^o)i and the Yukawa couplings h^: 

'i/'L means the SU{2) fermion doublet (^^)j and with Xr we denote the pair of 
the singlets ur, (Ir. 



2 AKY Configurations 

In order to minimize the energy of the static YM-H fields at a given value of 
the Chern-Simons number, 

^cs = -^J d'r E,,k {A^djAt + ^E^'^'^AfA'^Al^ , (2.1) 

one has to minimize the quantity 

£ = E,iass + 2^Ncs ■ (2.2) 
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Here C is a Lagrange multiplier and Edass is the energy of the static YM-H 
configuration which, in accordance with eq. ( |l.ip , is 

Eciass = J ci'rl^^iBtf + (A'J>)nA«>) + y (^^^'^ " t) } ' ^^'^^ 

where we have written Bf = ^^^ijkFjk- 

At the classical level the H field develops the v.e.v. (/) = {(()+) = 0, 

and the W boson field A'^ gets a mass mw = 9^/2, while the H field mass 
is rriH = Av. We choose the commonly used temporal gauge Aq[x) = and 
look for the minimum of eq. (|2.2| ) in a spherically-symmetric "hedgehog" form 



1,11: 



^"(r) = EaijUj h [dai - naUi) h Ualli 



^{v) = ^[H{r) + i{T.0)G{r)] |^ ° | , (2.4) 

where r = |r| and n = r/r. 

Gauge transformation which do not change the above ansatz are given by 

$ ^ C/$, Ai UAiU^ + iUdiU^ U = exp[i(n0)P(r)], (2.5) 

with Ai = AfT°'/2. Indeed, under this gauge transformation the functions A,B, 
C,H,G transform as follows: 

A{r) A{r)cos2P{r) - B{r) sm2P{r), 

B{r) B{r) cos 2P(r) + A{r) sin 2P(r), 

C{r) ^ C{r) + 2rP'{r), (2.6) 

H{r) H{r) cos P(r) - G{r) sin P(r), 

G{r) G{r) cos P{r) + H{r) sin P{r). 

This gauge freedom can be used to simplify the minimization of the functional 
( |2.2D for which we choose a gauge with C(r) = 0. This fixes the gauge trans- 
formation ( p.5D up to a global rotation with a constant P. 
We next introduce dimensionless quantities, 

gv 2 mir >. . 

x = r-=rm^, - = ^ = ^ = ^^'^^ 
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and will measure the energy in units of Mq = lirmw/a. The functional ( p.2| 
becomes: 



oo 

dx 



^2 ^ j^,2 ^ + B^^ 1)' ^ 2x2(G'2 + H'^) - 4BGH 

+2A{G^ - H^) + (1 + ^2 ^ B^){G^ + H^) + y ^^(G^ + - if 

+2CNcs ■ (2.8) 

The Chern-Simons number Ncs given by eq. ( p.l| ) is only well defined for boson 
fields which are continuous and vanish fast enough at infinity. In terms of the 
hedgehog ansatz ( |2.4| ) this means we have to use a gauge with ^(0) = A{co) = 
1, 5(0) = B{oo) = and C(0) = C(oo) = 0, which is not fulfilled for fields with 
C(r) = which minimize the functional ( p.SI ). Hence, before inserting ansatz 
( |2.4D into eq. ( |2.lD we have to gauge rotate the fields. A suitable transformation 
is obtained from eq. (|2.6[) with 5(0) = — 99(0)/2, P(oo) = —(f{oo)/2 where 
99(r) = arg [A(r) + ii?(r)]. By substitution of the gauge rotated fields into 
eq. (|2.1| ) we obtain in terms of the unrotated fields with vanishing C : 

Ncs = ^ J^^dx {A'B - B'A) + ^ (^^(oo) - (^(0)) . (2.9) 

Note that Nqs does not depend on the detailed shape of P{r). It is therefore 
possible to perform the minimization of £ with the unrotated fields A and B by 



inserting ( |2.9D into eq. (2^). Afterwards, the solution can easily be transformed 



into any other gauge by (2.6) 



In order to find the minimum of £ one has to solve the Euler-Lagrange eqs. 
for A,B,G,H with boundary conditions compatible with the finiteness of £. 
The details are given in app. A. One finds a family of solutions labelled by 
C; — 1 < C < 1. From eq. ( |2.9| ) one then gets the relation between Ncs ^-nd C 
which should be used to express E^ass through Ncs- The values of C = — 1, 
and 1 correspond to Ncs = 0, 1/2 and 1, respectively. 

If one is interested in the case of finite fermion number density, one has 
to add a term 2fiNcs to the energy functional ( |2.3| ), where fi is the chemical 
potential for fermions [Q, |5| . Since the form of the functional to be minimized 
is exactly the same as in the Lagrange multiplier approach, used to find the 
minimal energy barrier between the topologically distinct vacua, one proceeds 
in the same way as above. The only difference is that in the case of /x 7^ the 
energy is -E'dass(/^) = E^iassiO) + 2fj,Ncs ■ We present the minimal energy curves 
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as a function of the Chern-Simons number for different values of the chemical 



potential // in Fig. 1. At = we reproduce the curve of ref. [13|. It can be 
seen that at = Hcrit = ^ttuiw / a the energy barrier disappears - a fact which 
was anticipated previously from analyzing vacuum stability in the presence of 
finite baryon density |5[. At /x < ficrit the decay of a state with finite baryon 
density is given by tunneling "bounce" trajectories, found for small /i in ref. [||. 
It is a challenge for future work to determine the decay rate of a high-density 
state for any fi. 

One should not be irritated by the cusps of the curve at integer Nqs- the 
"coordinate" Ncs is a quadratic functional of the gauge field, and in a more 
natural coordinate the curve should be quadratic in the minima. To illustrate 
this, we have considered a 2-dimensional a model with a mass term, which 
mimics the 4-d electroweak theory. In this model the minimal-energy con- 
figurations can be found analytically (D.Diakonov, V.Petrov and M.Polyakov, 
unpublished), and the energy barrier is described by the simple but exact for- 
mula E = const IsimrNcsl, also exhibiting cusps at integer Ncs- However, 
if one sticks to just one degree of freedom (Nqs), one should ask what is the 
kinetic energy along this particular "coordinate". In the 2-d model this ques- 
tion can be answered exactly, and the restricted kinetic energy turns out to 
be Nl^sM{Ncs)/2 where the effective mass is M(Ncs) = const /\s'm7rNcs\, 
being singular at integer Ncs- Therefore, if one would replace Ncs by a coor- 
dinate X which would have a normal kinetic energy X'^ /2, the potential energy 
would behave normally at the minima ~ as X'^. 



3 Dirac Equation and Fermion Spectrum 

We are now starting to investigate the Dirac continuum and the fermion dis- 
crete level in the background field of the static minimal-energy configurations 
described in the previous section. 

We define the Dirac Hamiltonian as the operator TC appearing when one 



4-W)* = 0. (3.1) 



rewrites the Dirac equation following from eq. (1.2) in the form 

d 

Using the representation 

/=( : ° I , f =( " ? 1 , 75=1 1 , 0.2) 
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we can reduce the eight-spinors tp^, Xr to four-spinors tpi, Xr by the definition 

In the basis of the four-spinors ipi, the Hamiltonian TC is given by 

/ iaiDi M \ , ^ 

W= . . ' (3-4) 





where M is defined in eq. (|1.3| ). Here we have imphed the Aq = gauge. 

For the sake of numerical simphcity in finding the eigenvalues of Ti. we 
would like to exploit fully the spherical symmetry of the AKY configurations. 
However, if the masses of the upper components of the isospin doublets are 
different from the lower ones, i.e. hu ^ /id, the spherical symmetry of the Dirac 
equation would be spoiled. For that reason we shall only consider the case of 
equal masses, = = h. This is certainly a very good approximation for 
all leptons and quarks whose masses are much less than mw (which brings 
the scale into the problem) but it is not so good for the three (t, b) doublets. 
In Section 6 we will show how these massive doublets can be treated in the 
framework of this approximation. 

In the hedgehog ansatz ( ^.4| ) the matrices M, have the form: 

M = mp (H + iin0)G) hv , , 

, , mp = —p^ ■ (3.5) 

Mt = mF{H- i{n0)G) V2 

To find the fermion levels one has to solve the eigenvalue problem: 



(3.6) 



where 

^ ^ / icTiDi mp {H + i{n0)G) 

\ mp {H — i{ii0)G) —icidi 

It is easy to check that the Dirac Hamiltonian 7i ( p. 7] ) commutes with the 
so-called grand spin operator: 

K = L + S + t, (3.8) 
where L is the angular momentum, S is the spin, and T is the isospin operator: 

la ^ _^^abc^^Q^^ ga ^ l^a^ fa ^ l^a^ (3 9) 
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with = = 3/4. Therefore, TL can be diagonahzed in a basis of spherical 
harmonics with given grand spin K and its third component K^. According to 
the couphng rules of angular momenta, there are eight (in the case K = Q only 
four) basis vectors for fixed values of K, and radial momentum p (Both signs 
of the energy are included) . For the numerical diagonalization the basis is made 
finite by discretization of the radial momentum; the allowed values are denoted 
by pn- For details the reader is referred to app. B. We call the matrix elements 
of TL in this basis 'H'^Kj^,^{pm,Pn), where r, s = 1 . . . 8. Since they do not depend 
on we suppress this index in what follows; the degeneracy results in a factor 
2K + 1. The eigenvalues of TL'j^(pm,Pn) are denoted by exn or e^, the ones of 
the free Dirac operator H^^"^ are named e^^^ or e^^\ given by (e^)j)^ = + 

We have numerically calculated the eigenvalues function of Nqs 

for the transition from Nqs = to Nqs = 1- The most interesting behaviour 
is found in the K = sector, for which we plot some discretized eigenvalues 
in Fig. 2. At Ncs = (vacuum), there are no discrete levels and the spec- 
trum is symmetric; each level is two-fold degenerate. The levels of the lower 
continuum are occupied, the levels of the upper one are empty. When Ncs 
is increased, one level emerges from the lower continuum, travels all the way 
through the fermion mass gap, and finally reaches the upper continuum. This 
level stays occupied during the whole transition. As can be seen from Fig. 2., 
the other levels are slightly shifted upwards, each one replacing its predecessor. 
Eventually, at Ncs = 1 the spectrum looks exactly as at Ncs = 0, except that 
the first level of the upper continuum is now occupied, too. This shows that 
the transition between topologically distinct vacua is a baryon number violat- 
ing process. Moreover, for intermediate (non-integer) values of Ncs we have 
confirmed numerically the law B = Ncs following from the anomaly. Both 
the bound state level and the Dirac continuum contribute to the total baryon 
number B in this relation, and its proof requires an accurate gauge- invariant 
regularization of the Dirac sea, see Section 6. Previously, only the discrete level 
was investigated |16|; in this respect our results coincide with those of [p^]. 

Our next aim is to calculate the total energy contribution -E/er of the 
fermions. This energy is simply given by the sum of the energies of all oc- 
cupied states, relative to its value for the trivial vacuum at Ncs = 0. If we 
denote the energy of the discrete level which crosses zero by Evai, we obtain 



E 



Jer = E (^A - ) = E - E 4^^ + ^val 0{eval) , (3.10) 
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with the step function 6. Since Sp^H — Ti^^^j = 0, we can also write 

Efer = "9X^(1^^! ~ I^a'''|) + ^valf^i^val) = Egea + ^val (^{^val) ■ (3-11) 
^ A 

The sea energy Esea can also be written as 

Esea = -\ (spV?^^- SpVW) = -\Y.{2K + l)Y,(\eKn\ " \efS) ■ 

(3.12) 

The quantity -B/er is in fact a fermion one-loop correction to the classical energy 
of the boson field, in particular it gives the quantum correction to the sphaleron 
mass (another quantum correction arises from boson field fluctuations around 
the sphaleron). 

4 Renormalization 

The aggregate energy of the Dirac sea (|3.12| ) is divergent, since it is basically 
equivalent to the fermion loop in the external YM-H field. However, the elec- 
troweak theory is renormalizable which means that all divergencies can be ab- 
sorbed by local counterterms to the Lagrangian. There are exactly four types 
of divergencies corresponding to four possible counterterms: 

1. a divergency which can be absorbed in the counterterm, corresponding 
to the gauge coupling renormalization; 

2. a divergency which can be absorbed in the |D<l'p counterterm, corre- 
sponding to the Higgs "wave function" renormalization; 

3. a divergency which can be absorbed in the A^|<I>|^ counterterm, corre- 
sponding to the Higgs quartic coupling renormalization; 

4. a divergency which can be absorbed in the m|^|$p counterterm, corre- 
sponding to the Higgs mass renormalization. 

The last divergency is quadratic in contrast to the logarithmic divergencies in 
the first three cases. 

Any renormalization assumes a certain " regularization scheme"; we shall 
use a cut-off in a proper time representation for the Dirac sea energy. This 
quantity can be written as 

^sea = --SpV^=^| ^Spexp(-m2)^ (4.1) 
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where Sp is a functional trace, TL is the Dirac Hamiltonian (3.7), and r is 
the (inverse) ultra-violet cut-ofF, to be taken eventually to zero. A non-zero 



r suppresses Dirac levels with high momenta, \p\ > 1/\/t. In eq. (4.1) and 
below we do not write the subtraction of the free Dirac sea explicitly, though 
we assume it. 



Eq. (4J) implies that the classical energy functional (2^) contains bare 
couplings defined at the scale Xj^fr^ such as the gauge coupling (^(l/r), the 
Higgs expectation value T;(l/r), etc. However, when one finds e.g. the classical 
mass of the sphaleron, one expresses it in terms of the physical couplings and 
Higgs expectation value, defined at the electroweak scale like mvi/- Therefore, 
in order to use the standard values for the physical constants, a proper time 
integral between r and Xjm^ has to be absorbed by the classical energy. To 
be more precise, we write: 



dt 



Sp exp(— i'H^ 



+ 



(4.2) 



In the first integral one can perform a semi-classical expansion as it deals 
with the spectral density of the operator H at high momenta. The technique 
is explained e.g. in ref. [|l^]. We get for the divergent piece of the Dirac sea 
energy: 



1 



4v^ 



^ Sp exp(-m^) 



1 



167r2 



dt \ d-^r ■ 



2g m 



m 



^(cDtcD-IL) 



w 



1 



1 



4V^ 



'"w 

'^w dt 



+ 0(t° 



Sp exp(— t'H^ 



div 



+ finite terms. 



(4.3) 



We see that the quadratically (1/t^) and logarithmically (1/t) divergent 
terms are exactly those entering the classical energy functional ( |2.3| ). Therefore, 
they can and should be combined with the bare constants of the corresponding 
terms in the classical energy - to produce renormalized physical constants at 
the scale of my/. For example, the coefficients in front of the magnetic field 
energy, -B^/2, combine into 

1 1 , . o . 1 



(4.4) 



52(l/r) ^SiT-^-^-^'' g^irnl,)' 

and similarly for the other three terms. (Eq. ( |4.4| ) is not the whole renormal- 
ization of the gauge constant since boson fiuctuations also contribute to it). 
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Because the divergent part of ( |4.3| ) is combined with the classical energy, the 
renormalized Dirac sea energy is: 

KTa = ^ Sp exp(-tW2) + y^^^ ^ Sp exp(-m^) 

'"^^ '^^ [Sp exp(-m2)l \ . (4.5) 



div 



This expression is finite as r ^ 0, and one can put r = in it. Introducing a 
dimensionless variable tm^r — > t and measuring all quantities including the sea 
energy in units of mw we get: 

Though eq. ( [4.6| ) solves the problem of renormalization it is still not too 
practical from the numerical point of view. The divergence of the unrenormal- 
ized Dirac sea energy is due to the large contributions for high values of K. 
Therefore it is convenient to perform the renormalization subtraction for each 
K sector separately, like 

Eli: = Ei^K + 1){Ek - Ef") . (4.7) 

K 

Then at each step one deals with perfectly finite quantities, and the summation 
in K is convergent. 

To this end we calculate the divergent part ( [4.3| ) of the functional trace (i.e 
the quantity to be subtracted) in a basis of a complete set of eigenfunctions of 
the free Dirac Hamiltonian, given in app. B. From the definition of the trace 
we have: 



K 

rR 



•V V / dxx'^vP-{pn,x) exp(-ml^) ..u^\(Pn,a;) . (4.8) 

r=l 71=1 

where TCk is the Hamiltonian in the basis of the spherical eigenfunctions for 
given grand spin K and radial momentum p„, in fact an 8 x 8 matrix, and 
v^^\{pn,x) are the radial parts of the eigenfunctions. For details see app. B. 

The semiclassical expansion in a given K sector corresponds to expanding 
the exponent in eq. (4.8) in powers of TIk —Pn- The expansion is rather tedious 
and the result is lengthy. We present the final formula in app. D. There we also 
show that eq. ([4.3[) is reproduced from eq. (WM) in the limit R ^ oo. 
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5 Renormalization at Non-Zero Tempera- 
tures^ 



Keeping in mind the application to the baryon number changing sphaleron 
transitions at non-zero temperatures (below the electroweak phase transition) 
we explain here how the renormalization is performed for the baryon number 
changing rate. This question has been addressed in ||l^ and [^] for the limit 
of high temperatures; it was demonstrated there how the theory is reduced 
to three dimensions. In our work we will not assume the high temperature 
limit but consider temperatures in the range ( |5.2D (see below). This leads to 
additional contributions which are not present in the limit of high temperature. 

In principle we treat both fermion and boson fluctuations on equal footing. 
Our numerical calculations, however, are restricted to the corrections due to 
fermion fluctuations. The bosonic part will be numerically dealt with in a sub- 
sequent paper; so far it has been investigated in e.g. |^]. In an abelian 
(1 + 1) dimensional model it is possible to calculate the sphaleron transition 



rate analytically. This has been done in |23] for boson loop corrections and in 
[ p^ ] for fermions. 

The traditional starting point is the semi-classical Langer-Affleck formula 



26, ETl for the thermal transition rate F: 



r = ^ImF, (5.1) 

TT 

where a;_ is the negative boson eigenfrequency around the sphaleron, /3 is the 
inverse temperature and F is the free energy of the system]^ One can imagine 
introducing a small chemical potential for baryons so that the minima at integer 
Ncs become metastable, and \mF gets a precise meaning Q. 

The Langer-Affleck formula works in the temperature range given approx- 



imately by 



'^<T<E:,Zs{T). (5.2) 



^This section was prepared in collaboration with V.Petrov and P.Pobylitsa. 

■^A question may arise about the applicability of this formula to the process under consideration 
since it implies that the transition occurs mainly above the sphaleron barrier. It is not too clear 
whether it is justifiable to treat such degrees of freedom as belonging to thermal equilibrium, and 
whether presumably perturbative field fluctuations above the barrier do actually lead to the baryon 
number change. Keeping in mind these reservations we nevertheless start from eq. ( ^.l[ ). 
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If the temperature is lower than u;_/27r, the baryon number violation rate 
is dominated by transitions via periodic solutions of the classical Yang-Mills 



equations in Euclidean space [28|. In the case T = these solutions reduce 
to instantons. Numerically the value of u)-/2tt was found to be « 0.3 miy for 



niH = mw [20, 29 1 . The upper bound E'^[assi'^) renormalized energy 

of the sphaleron barrier. As we show below, it depends on T; it vanishes at 
some critical temperature Tc where the symmetry breaking of the electroweak 
theory disappears. Hence, eq. ( ^I^ ) implies T < Tc. We will discuss numerical 
results in Section 6, but let us mention here already that for mn = mw and a 
reasonable value of the top quark mass, Tc is of the order my\r to 1.5 mvK- We 
see thus that the domain of applicability of eq. ( |5.lD is rather restricted, and 



the reduction of the problem to three dimensions |20|, |21|, ^ implying T » mw 
does not seem to be justified in reality. 

In the semi-classical approximation one takes the static sphaleron configu- 
ration as a saddle point for ImF, and has to compute the boson and fermion 
determinants around it. Special care must be taken towards the zero and nega- 
tive modes of the boson small-oscillation operator; this problem has been solved 



in ref. |2C, 21], and we denote the resulting factor by A/q,-. The boson deter- 
minant below involves only positive eigenmodes, and is denoted by a prime. 
Eq. (p.lp can be written as 



where "bare" means that the classical mass of the sphaleron is computed with 
bare physical constants defined at some ultra-violet cut-off scale. The super- 
script (0) refers to the free (no field) determinants. 

Let us denote by the eigenvalues of the 3-dimensional quadratic form for 
boson fiuctuations about the sphaleron and by e„ the eigenvalues of the Dirac 
Hamiltonian for fermions; the same quantities with a superscript stand for 
the corresponding no-field eigenvalues. Imposing boundary conditions periodic 
in time for bosons and antiperiodic for fermions, one immediately establishes 
the eigenvalues of the 4-dimensional small-oscillation operators {k G 2): 

ASf=HU(^)'. A/;=.„+,:fi^. (5.4) 
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The determinants of eq. (^.3D are the products of these eigenvalues: 



Det 



bos 



A 



bosfl 
n,k 



,Det 



(0) 



bos ' 



n n ^=n 



Det 



fer 



Det 



(0) 
fer' 



n n 



\fer 



n k=—oo ^n,k 



ferfi 



n 



sinh(a.O/3/2) 
sinh(a;„/3/2) 

cosh(e„^/2) 



cosh(e0^/2) ■ 



(5.5) 



At this point McLerran et al. M, 20, 21 1 argue that in the case T ^ rny^/, one 



should replace sinh by its argument and cosh by unity, since the 3-d eigenvalues 
Wri, £n are of the order of mw- We will show later that this replacement is 
justified, but as the condition T ^ mw is not compatible with the allowed range 
( |5.2D we have to use the exact expressions in what follows. The determinants 
(|5.5D can be identically rewritten as: 



sinh(a;0/3/2) 



n 



sinh(a'„/?/2) 

cosh(e„/3/2) 
cosh(eO/3/2) 



exp 



p 




+ Eln 



(5.6) 



(5.7) 



The first terms in the exponents of eqs. {^^, 5.7) make up the zero-point oscil- 
lation corrections to the bare sphaleron mass at zero temperatures: 



^cra7. + E^(^n-'^°)-^E(l^ 



2' " 2 

n n 



E^ 



(5.8) 



Both bosonic and fermionic sums are divergent but the divergencies are ab- 
sorbed in the renormalization of the constants entering E^f^^^ so that Ei is 
finite. In Section 4 this renormalization was explicitly performed for the case 
of fermions. 



The terms containing logarithms in eqs. (5^, 5^) depend on the tempera- 
ture but are ultra-violet finite, since the large eigenvalues are cut by the Boltz- 
mann factors. For a closer inspection of these terms we use spectral densities, 
in which the subtraction of the vacuum is included: 



1 



) - ln(l - e 



1 , 

dE p^"'{E) In 



1 



-I3E 



(5.9) 



E [ln(l + e-^l^"l) -ln(l + e-^l^"l)j = H dE p^^^ (E) In (l + e'l^^^ 

„ -' J— oo ^ 



(5.10) 
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We can derive a semiclassical high-energy expansion for the density, which 
is related to the expansion of eq. If the latter one is written as 

Sp (^e-*«' - e-*«'°'') = at-V2 ^ ^^1/2 ^ ^^3/2 ^ ^ (g ^^) 



one finds: 

pf^'iE) = \E\ Sp (<5(W2 _ _ 5(7^(0)2 _ ^2 

1 / 5 1 3c 1 \ , , 

We see that the first term corresponds to a constant spectral density p^** at 
large E (which, of course, is in perfect correspondence with the quadratic di- 
vergency of the Dirac energy). We get from eq. (4.3) 



and a similar expression (but with a different coefficient) for p^'^. Putting these 
Poo into eqs. ( p.9| , 5.101 ) we obtain for fermions: 



(5.14) 

This expression is obtained for each fermion doublet so that we have to sum 
over all doublets. Additionally, we must consider a similar result for bosons, 
where mjp has to be replaced by (3/4)(m^ + 3m^) in the prefactor |3^. At 
T = 0{m\Y/g) these terms can be of the same order as the (renormalized) 
classical zero-temperature energy of the sphaleron (divided by T), therefore in 
that range of temperatures one has to find a new sphaleron solution, with a new 
functional which includes the thermal part ( |5.14| ) . Fortunately, it has the same 
form as already encountered in the classical energy functional (see eq. (p.3|)); 
therefore, instead of solving anew the classical eqs. of motion one can use the 
non-thermal solution for the YM-H fields but with the parameter v replaced 
by the temperature dependent expression 

^(2^)2 ^ ^2 _ 3m^ + 9m|t, + 4Edoubl. ^2 _ ^^^^^^ 

This leads to a temperature dependent renormalization of the masses: 

mniT) ^ mw{T) ^ mpiT) ^ v{T) 
niH mnw "nip v 
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With this thermal renormahzation of the parameters performed, one has to 
subtract the quantities p^^'^^''"' from pbos, fermi^^^ g^g^ ^ ). After 
the subtraction high E are suppressed as in the integrand so that in the 

hmit T ^ niw the contributions from eqs. (|5.9| , 5.10I ) read: 



dE[p^'^^{E)-p''^)\u{m. (5.17) 

r di?(p/-(i?)-/,Mln2. (5.18) 

Apart from the subtleties with the renormahzation mentioned above, this 
is what one would naively get from eq. ( |5.5| ) by replacing the fermionic cosh by 
unity and the bosonic sinh by its argument, which corresponds to the recipe 
of ||6|, ^0|, However, we are not working in the high temperature limit but 
in the range eq. ( |5.2|) so that we have to use the full formula for the sphaleron 
transition rate: 

r = I^ATo,- exp(-/3^tot(r)) , (5.19) 



where 



with 



E,^m = E^ZsiT) + C7 + > (5-20) 

POO , \ F 

KZ" = dE (p^-(i?) - - 0{1/E')) I 

+ HdE [p'"^%E) - p'Z^) ln(l-e- 

/oo . . 

_JE{pf^^iE)-pt-0{l/E')] 

- dE (/-(i?) - pt) In (l + e-'^l^l) . (5.21) 



\E\ 



Here E^[^^g{T) is the classical sphaleron mass with the temperature-re- 
normalized values of the physical constants given by eqs. ( ^.15 , ^.16 ). According 



to eqs. ( 5.20 , f).2l\ ) the classical mass gets corrections due to the temperature 



independent zero-point oscillations (the first terms of E^^^^ and E^j^^) and 
due to temperature dependent terms (the last terms of E*^^^ and E^^^^). The 
0{1/E'^) terms to be subtracted from the spectral densities correspond to the 
logarithmic divergencies written explicitly for the fermions in eq. ( 5.12| ). 



The boson loop corrections I3E^^^^ are suppressed by the factor a so that 
they are usually discussed in context with the pre-exponential factor. Its in- 
vestigation is not performed in this work. The fermionic contribution (3E^j^^ 
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is also suppressed by a, but it is enhanced by the number of fermionic species 
Nf = 12 so that its value can be considerable. We shall see in next section that 
it suppresses the transition rate significantly, especially for massive fermions. 



6 Numerical Results and Discussion 

In this section we present our numerical results. Since the dependence on the 
Higgs mass rriH turns out to be weak we decided to fix mn = rny^/ for all 
calculations presented in this section. We will briefly comment on the influence 
of ruH at the end of the section. All results are computed using the physical 
value g = 0.67 for the coupling constant. 

Let us start by discussing the behaviour of the fermion spectrum and the 
contributions to the energy for the path from Ncs = to Ncs = 1- The level 
crossing phenomenon has already been demonstrated in Fig. 2. Additionally, 
we have investigated the discrete level for various fermion masses mp and our 
results agree with those of Kunz and Brihaye |l^ . 

In the following we present numerical results for the fermionic energy -E/en 
especially its sea part -EJI" , given by eqs. ( p.lO| , p.ll| ) and eq. (^). Eq. 



implies that the UV cut-off A = \l \/t has to be taken to infinity. In practical 
terms, however, we always have to work with a finite A since our numerical basis 
must be finite. Therefore we calculate i^^ea ^-s a function of A and extrapolate 
to its value £'^|^(oo) at A = oo. From the semiclassical expansion of i?sea for 
high A we obtain the following dependence: 

EITM = (oo) - ^ + O (^-1) , (6.1) 

with some constant C . (The terms ~ A^ and ^ ln(A) are removed by renor- 
malization.) 

Apart from working with finite A it is also necessary to introduce two nu- 
merical parameters, the box size R and the maximum momentum Pmax (see 
app. B), in order to obtain a finite basis for the diagonalization of H. For a 
fixed value of A both have to be taken large enough so that E^^^i^K) does not 
change any more when they are further increased. In Tab. || we plot results for 
fixed rap = 2mwi ^cs = 0.5 and A = 5mw for different R and Pmax- It can 
be seen that stability is reached for R ~ 12 and Pmax ~ 3A so that we 
obtain as our result: 

Ell2i^ = ^rnw,mF = 2mw,Ncs = 0-5) = 8.55 mw, (6.2) 
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P /A 


1.5 


2.0 


2.5 


3.0 


3.5 


Rmy/ = 5 


7.98 


7.29 


7.20 


7.20 


7.20 


Rmy/ = 7 


8.99 


8.34 


8.25 


8.25 


8.25 


Rmw = 9 


9.25 


8.57 


8.49 


8.48 


8.48 


Rmw = 10 


9.25 


8.61 


8.52 


8.52 


8.52 


Rmy/ = 11 


9.28 


8.62 


8.54 


8.54 


8.54 


Rmy/ = 12 


9.31 


8.64 


8.55 


8.55 


8.55 


Rmy/ = 13 


9.31 


8.64 


8.55 


8.55 


8.55 


Rmy/ = 14 


9.28 


8.64 


8.55 


8.55 


8.55 



Table 1: £'^|"(A = bmy/) for different values of the numerical param- 
eters R and Pmax with mp = 2mw, Ncs = 0.5 and A = 5mw fixed. 
Stability is reached at 12m^ and Pmax ~ 3 A. 

with a numerical error of less than 0.2%. Using the same method we obtain 
results for other values of A and for different mp and Ncs- It turns out that 
in general it is sufficient to choose R = 12 and Pmax = 3.5 A. In Tab. ^ 
we present results for various values of A with R = 12 m^ fixed. Additionally 
we have written values for £^J|^(A) obtained from eq. ( |6.1D with El^'^(oo) = 
9.09 my^ and C = 13.54 m|^^ into the last column. A comparison with the 



Pmax / ^ 


1.5 


2.0 


2.5 


3.0 


3.5 


from eq. (0|) 


K/my/ = 3 


8.28 


7.71 


7.64 


7.63 


7.63 


7.59 


A/mvF = 4 


8.86 


8.33 


8.25 


8.25 


8.25 


8.24 


k/my/ = 5 


9.27 


8.63 


8.55 


8.55 


8.55 


8.55 


k/my/ = 6 


9.54 


8.80 


8.71 


8.71 


8.71 


8.71 


k/my/ = 8 


9.89 


8.98 


8.88 


8.88 


8.88 


8.88 



Table 2: i?^|"(A) for various values of A and Pmax- R = 12m^\ 
mp = 2my/ and Ncs = 0.5 are fixed. The values in the last column 
are obtained by eq. (6.1) with EH^^oo) = 9.09 m^y and C = 13.54 m^. 

numerical results shows that the power law eq. ( |6.1[ ) is fulfilled very accurately. 
Therefore the numerical error of the extrapolated value -E'g|"(oo) is of the same 
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order as the error of ii^J|"(A) for finite A. Thus our extrapolated result for -EJI" 
finally reads: 

El^'^^imp = 2mw,Ncs = 0.5) = 9.09 mw , (6.3) 

again with a numerical error of 0.2 %. Our results for different values of and 
Ncs were obtained the same way. The numerical errors always have roughly 
the same size. The relative error can increase up to 2% for small mp where -EJI" 
is considerably lower than for nip = 2mw- This accuracy is fully sufficient for 
our purpose; one should keep in mind that the error arising from taking equal 
masses for the top and bottom quark is surely bigger than the numerical error. 

In Fig. 3 we have plotted E'p,^ = -E^l" + s^ai d{£vai) as a function of Ncs 
for mp/mw = 0, 1 and 2. For mp = there is no contribution from the 
valence level, the energy is symmetric with respect to Ncs = 0.5 . Its total 
value is very small (-E^la = 0.18 mi^/ for Ncs = 0.5) compared to the classical 
energy of the boson fields (ps 100 m^y for Ncs = 0.5). This is also the case for 
mp = m\Y , where the energy is the sum of the symmetric sea and antisymmetric 
valence part. The fermionic energy Eyj^ is only significant in comparison to the 
classical energy when mp \s much larger than mw ■ For mp = 2 mw the energy 
of the fermions is already about 9% of the classical energy, and it increases 
dramatically with mp. In fact, one can show that in the leading order the sea 
energy grows with mp as mp ln{mp). This law is already a good approximation 
for mp > 3m\Y- In Tab. ^ we have plotted results of -E^ea different mp at 
Ncs = 0.5 . 



mp/mw 





1 


2 


3 


4 


5 


6 


E7ea/mw 


0.18 


0.51 


9.05 


72.63 


303.5 


862 


2001 



Table 3: -E^l" as a function of mp for fixed Ncs = 0.5 . The energy 
increases according to ~ mpln{mp) . It should be compared to the 
classical energy which is about 100 for Ncs = 0.5 . 



Before we turn to the discussion of the temperature dependent contribution 
E^j^^ to the sphaleron transion rate and the baryon number and -density, let us 
briefly comment on the numerical accuracy of these UV- finite quantities. The 
only parameters which turn up in the numerical evaluation are R and Pmax- 
Again, both have to be chosen large enough to ensure stability of the results. 
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This is done in analogy to the method for the calculation of £'J|"(A) for fixed 
A. The numerical errors are again in the range of 1% to 2%. 

We are now going to present the temperature dependent contributions to 
the sphaleron transition rate, given by eqs. ( P1^5^ ). In this paper we only 
considered the classical term /^-Ecfoss ^'^d the fermionic part f3E^j^^ but not the 
bosonic part /JE"^^™^. The temperature enters the transition rate in three ways: 
First, in form of the prefactor (5 = 1/T which causes the transition rate to 
become unsuppressed when EtotiT) T. Second, the masses get renormalized 
via eq. ( 5.16 ) so that the energy EtotiT) is roughly proportional to v{T)/v(fi) 
given by eq. (|5.15| ). Third, the last term of E^^^^ is explicitly dependent on (3. 
The first two effects are numerically trivial but very large, while the last one is 
difficult to compute and rather small. 

We start by considering the first two aspects. The dominant contribution 
in eq. ( |5.2[1| ) is 

For r ^ this terms goes to infinity, but it decreases rapidly with increasing 
T. In the vicinity of the critical temperature 

'"'"^^ (6.5) 



?,m]j + 9m^ + 4 Edoubl. m'p ' 

where the symmetry breaking of electroweak theory and the sphaleron barrier 
disappear the function /3£'^fJ^^ (T) goes to zero and the baryon number violation 
rate becomes unsuppressed. 

The value of Tc and the quantitative behaviour of PE^f^^^ (T) depend largely 
on the choice of the parameters. We wish to focus our interest to a situation 
which resembles the physical one as close as possible. As it was mentioned 
above to preserve spherical symmetry of the Dirac equation one has to consider 
the case of equal masses for up and down fermions. This approximation is 
not so good for the three (i, b) doublets. In the case rrit^b < it^w the fermionic 
corrections can be estimated using the average mass of the doublet, and treating 
the difference as a perturbation. The corrections are anyhow small in this case. 
In the more realistic case mf, <^ m\Y <C mt the correction from the top quark 
is half that of the doublet with both masses equal to nit (This recipe can easily 
be derived from eqs. ( |4.5| ) and ( [5. 12]) in the limit nit ^ mw)- Therefore, in our 
numerical estimates we take 9 + 3/2 massless fermion doublets and 3/2 massive 
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doublets with mass mt which we vary in the range rrit/mw = 1-5, 2.0, 2.5. In 
particular the sums in eqs. ( |6^ , |6.5| ) are replaced by 3/2 m^. 

We have plotted PE^iassC^) ^ ^ function of T in Fig. 4 (dashed lines) for the 
masses mt/m\\r = 1.5, 2.0 and 2.5. We realize that Tc is rather low (between 
mw and \.f>m\\r). One can obtain a higher by taking a larger Higgs mass, 
but this will not change the picture qualitatively. The current experimental 
bounds do not suggest that mn should be considerably higher than mw- We 
have already stated in Section 5 that the range of applicability of the Langer- 
Affleck formula is roughly given by 0.3 m^^ ~ u;_/27r < T < Tc ~ l.f>m\\r. This 
is the range we use in Fig. 4, but one has to be careful if T is very close to Tc- 
In this case the Langer- Affleck formula is clearly not valid because the barrier 
gets so low that the transition occurs far above the sphaleron. 

We now want to investigate how far the fermionic fluctuations (3E^j^^ of 
eq. ( ^.21 ) influence the transition rate. The numerical results for this contribu- 
tion are included in Fig. 4 (solid lines), again for the three choices rrit/mw = 1.5, 
2.0 and 2.5. We find that the results for these cases are quite different. For 
mt = 2.5 mw and T ~ 0.5 m^y the fermionic fluctuations increase the exponent 
of the classical Boltzmann factor by about 50%. For T ~ 0.9 mw this increase 
is even about 65%. On the other hand, for mt = 1.5mw the fluctuations are 
only between 5 and 9% of the classical exponent. We conclude that the fermion 
loops suppress the baryon number violation rate in any case, but the strength 
of the suppression depends on the value for the mass of the top quark. We 
expect that the above numbers slightly change if the difference between mt and 
mb is treated perturbatively in higher orders. 

The suppression of the transition rate is important for the understanding 
of the matter excess in the universe. However, such cosmological consequences 
are still a matter of discussion (see e.g. |lO[), and we will not go into details 
here. 

Let us turn to the discussion of the baryon number and -density. First 
we have to find reasonable definitions for these quantities. In principle, the 
baryon number is the number of all negative energy eigenstates of the Dirac 
Hamiltonian, relative to the vacuum. The valence state always contributes to 
the baryon number, so it has to be added explicitly after it has crossed zero. 
However, the increase of the baryon number according to the law B = Ncs is 
not only caused by the level crossing of the valence state, but it is also due to 
a process that happens at both ends of the spectrum. To see this let us return 
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to Fig. 2 which demonstrates that, going from Ncs = to Ncs = one state 
emerges from "outside the spectrum" and joins the negative continuum, and 
another one disappears "from the positive continuum to infinity". In order to 
take this effect into account correctly we have to introduce a regularization 
and compute the result with a finite cut-off which eventually must be sent to 
infinity. Thus we define: 

B = lim f N{ex < 0) - A^(ef < 0)| + 9{e,ai) 

A^oo I- "J reg 

= lim < 

A^oo 2 




= lim -5^-^ Tdtt-i/^e-^A* (6.6) 

It is possible to show analytically that 5B = 6Ncs with Nqs and B being 
defined by eq. (|2.l| ) and eq. (|6.6| ), respectively. Here "5" means the variation 
with respect to the boson fields. One has to differentiate the integral in eq. ( |6.6D 
with respect to e\, perform the t integral and expand in powers of 1/A. The 
only term which does not vanish in the limit A ^ oo is identical to 6Ncs- For 
the baryon density we obtain: 



+ 9{e,ai) f dnipl^i{r)ip,ai{r). (6.7) 



Here if\{r) are the eigenfunctions of the Dirac Hamiltonian: Tiifx = e\Lp\. Of 
course, by definition we have B = drr'^p{r). 

In order to evaluate eqs. ( |6.6| , |6.7| ) numerically we have to insert the eigen- 
values Ex obtained from the diagonalization in different sectors of grand spin 
K. It turns out that B and p are dominated by the K = sector. In fact, its 
contribution is usually more than 90% of the total value. This means that the 
created baryon is basically in a = state. 

In Tab. ^ we present the baryon number B as a function of Ncs for mp = 
mw- We observe that the law B = Ncs is excellently reproduced. The same 
behaviour is obtained for different masses mp. The baryon density as a function 
of the radial distance r is plotted in Fig. 5 for Ncs = 0.26, 0.5, 0.74 and close 
to 1 {niF = mw). Between Ncs = and Ncs = 0.5 the density increases 
with Ncs for each fixed r. For Ncs > 0.5 the density at the origin p{0) 
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Ncs 


0.00 


0.15 


0.26 


0.35 


0.43 


0.50 


0.57 


0.65 


0.74 


0.85 


1.00 


B 


0.00 


0.18 


0.27 


0.35 


0.43 


0.50 


0.57 


0.65 


0.73 


0.82 


1.00 



Table A: B as a function of Ncs for nip = niw- The data confirm the 
law B = Ncs coming from the anomaly. 



decreases again, while the integral B is still increasing (as seen from Tab. 
Thus we are creating a baryon which becomes more and more delocalized. 
Approaching Ncs = 1 the density becomes independent of r and infinitesimally 
small everywhere. This corresponds to a free, unbound particle as expected 
from a vacuum configuration of the classical fields. 

Finally we have checked the influence of the Higgs mass mn on the preceding 
results. Qualitatively, the picture was found to be the same for any mn (at 
very large mn, however, the sphaleron solution becomes modified, see [p!^). 
In fact, we observe that both and E^j^'^ are slightly decreasing functions 
of mn- So for mu > niw they always have the same order of magnitude 
as their classical counterparts. Only for the unphysical value mn — > the 
fermionic terms show a significant increase. In this case the energy i^jg" reaches 
approximately 10% of the classical energy {mp = myy), which is about 20 
times higher than for mu = m\y. This relatively large value is due to the slow 
decrease of the boson fields for r — > oo which causes a strong perturbation of 
the Dirac sea. 

Acknowledgement: We are grateful to V.Petrov and P.Pobylitsa for nu- 
merous discussions. D.D. and M.P. would like to thank the Institute for Theo- 
retical Physics II of the Ruhr Universitat Bochum for hospitality. Their work 
was sponsored in part by the Russian Foundation for Fundamental Research 
under grant 7,^:93-02-3858 and by the Deutsche Forschungsgemeinschaft. 



Appendix A 



For the sake of completeness we cite the eqs. for getting the minimal-energy 
configurations of Akiba et al. |l^]. They are found from the Euler-Lagrange 



eqs. which follow from varying eq. (2^ 



A" = \a{A^ + B^-1)+A{H^ + G^) + G^-H^-2Cb', 



24 



B" 
{xHf 



2x 



;B{A^ + 52 - 1) + 5(^2 + G2^) _ + 2Cy4', 



H[{A-lf + 3"^ 



\bG+^xH{H^ + G^-1), 
X I 

.2 



(xG)" = —G[{A-lf+B^] + -AG--BH + —xG{H'^ + G^-l) 
2x X X 2 



(A.l) 



A consequence of these eqs. is a "conservation law" which can be used to check 
the accuracy of the numerical integration: 



{AB' - A'B) + x'^iHG' - H'G) + C(l - A^ - B^) 



0. 



(A.2) 



The boundary conditions for eqs. (A.l) follow from the requirement that 
there are no singularities at the origin and that the energy functional does not 
diverge at infinity. The first requirement implies the following behaviour of the 
functions near the origin: 



A 
B 
H 
G 



cos a (1 + rx'^) — |(st + 2rC)(sin a) x^ + ... , 
sin a (1 + rx^) + i(st + 2rC)(cos a) + ... , 
s cos Y + t{sm f ) X + j2n'^s{s'^ — l)(cos x^ + 
s sin ^ — t(cos f ) X + |^K^s(s^ — l)(sin ^) x^ + . 



(A.3) 



where a, r, s, t are arbitrary constants. 

Requiring the convergence of the energy integral at infinity and solving the 
eqs. (A.l) at large x one gets the following behaviour of the functions at x ^ 00: 



A 
B 
H 

G 



cos 7 + e sin(Cx + /3) -E'(x) 
sin 7 — ecos(Cx + /3) -E'(x) 



h 

cos TT I 1 H — e 

X 

sm ^ ( 1 + — e~ 



' j - e sin ^ cos(Cx + /? - 7 + (5) 



+ e cos ^ cos 



(Cx + /3-7 + 5)^, (A.4) 
x^ 



where we have denoted 



E{x) = exp (^-\/l-C^x^ , C 



sm 



2 ' 



1 — = cos I , 



and 7, P, e and h are arbitrary constants. 



(A.5) 



Note that the asymptotics ( A.4 ) is valid only if rriH ~ mw, if the Higgs is 
much heavier than the W boson, eqs. ( A.4 ) have to be modified. 
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We have thus four arbitrary constants defining the behaviour of the func- 
tions at infinity 7, e, h) and four in the origin (a, r, s, t). Matching solutions 



of eqs. (A..1) starting from the origin with those starting from infinity, one finds 
these constants for a given value of C (and k) - up to a global gauge rotation, 
however. The point is, only the combination 7 — a is invariant under the gauge 
rotation (|2.6D with a constant P. One can use this freedom to fix a = 0, for 
example. 

It should be noticed that eqs. ( |A.1| ) have a solution only for |C| < 1. The end 
points of this interval correspond to Ncs = 0, 1; it is seen from eqs. (A. 4, A. 5) 
that at |C| — > 1 the functions rather oscillate than decrease. 



Appendix B 



In this appendix we summarize how the Dirac Hamiltonian 7i of eq. (3.7) is 
diagonalized. Since it commutes with the grand spin operator K, we diagonalize 
it in a basis of eigenfunctions of and K3. To this end we introduce normalized 
eigenstates of the commuting operators L^, = (L + S)^,K^ = (J + T)^ and 
its third component K3. Naturally, they are also eigenstates of the operators 
g2 _ -ji2 _ 2^^^ shall use the short-hand notation \K, J, L; K-^) to denote 
these states. 

Since only the operators K^, K3 commute with the Hamiltonian Ti but 
not or L^, states labelled by different J,L at given K, will mix in the 
eigenvalue eq. (^^). Generally speaking, for given values of K, the following 
four states mix: 

\K, K + ^, K+1;K3), \K, K + ^, K;K3), 

\K, K K;Ks), \K, K K - I; K3). (B.l) 

We therefore have to decompose the spinor-isospinor functions i/jlj^r in the 
above four states for every value of K, K3: 

00 K J+i 
K=QK^=-K j=K-\L=J-\ 

00 K J+h 
^R(r) = E E E E ^''9'^^'{r){m,J,L■,K,). (B.2) 
A'=0 K-A=-K .J=K~^ L=J-^ 
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The eigenvalue eq. ( |3.6| ) is thus mixing four functions / (from ipi) and four 
functions g (from Xpt) for any value of K, K^. Since the Hamiltonian is block- 
diagonal with respect to K^,, eq. ( |3.6| ) results in a system of eight ordinary 
differential eqs. for eight independent functions for each K. In the case of 
K = {) where the last two states in eq. ( p3.1D are absent, only two states for ipL 
and two for X/j mix, and the eigenvalue equation involves only four functions. 
This particular case was considered recently in ref. 

The set of differential equations will be solved by putting the system into a 
large but finite spherical box. To insure that at the boundaries of the box the 
boson fields take their vacuum values, we have to switch to a gauge in which 
A{oo) = 1, B{oo) = 0, C(oo) = 0, H{oo) = 1 and G(oo) = 0. Moreover, to 
insure continuity at the origin the values of the fields at r = must be A(0) = 1, 
-6(0) = 0, C(0) = and G(0) = 0. Therefore we have to give up the simplifying 
choice C(r) = and to include the C-field into the formulae. This is, however, 
no substantial difficulty. 

In substituting the decomposition ( p. 21 ) into eq. ( |3.6| ) one has to calculate 
matrix elements in the \K^ J, L; K^) basis of the scalar operators appearing in 
eq. (p. 61), namely (osn), (o90), (0n), (cef?), ([cex0]n) and (oen)(0n). We call these 
operators "scalar" as they commute with the operator K^, hence equations for 
different K will not be mixed. 

Thus, the Dirac Hamiltonian ( |3.7D is block-diagonal in the \K^ J, L; K^) ba- 
sis, where the blocks correspond to sectors with definite K and X3; we call these 
blocks TCkKs- Since they do not depend on this index will be suppressed 
subsequently. Each Tix is a 8 x 8 matrix acting on an 8- vector Vk- 

HkVk = EVk, (B.3) 

where 

9k+^,k+v 9k+^,k^ 9k-\,k^ 9k-1,k-i)^ (B-4) 
and Ti.K can be presented as 

[ -Vk 
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Here Vk is a 4 x 4 matrix composed of differentiation operators: 



\ 



_d_ _|_ K+2 

dx ' X 






dx X 











A. -I- K+i 

dx X 







_d_ I K-l 

dx X 



\ 







(B.6) 



while Vk and Wk are also 4x4 matrices made of the background W and H 
fields. To shorten notations we introduce 



Ak 



1 - A{x) 

XCK 



Bk- 



K{K + l) 

B{x) 



XCK 



CK -- 

C{x) 

XCK 



2K + 1, 



Gjx) 

CK 



Ck = \ ' , Gk = ^ ' , Hk = H{x). 

(B.7) 

Using standard 6j symbols as described in app. C, we obtain for the matrix 
Vk: 



I Bk{K + 1 

-Ak{K + 1) 

AKbK 
\ {Bk - CK)bK 



Ck 
2 



Ck 
2 



and 



Wk = nip 



-Ak{K + 1 
-Bk{K + 1)- 
{Bk + CK)bK 
-AKbK 



Hk Gk 

—Gk Hk 

2GKbK 

2GKbK 



AKbK 

{Bk + CK)bK 
-BkK + % 
AkK 

-2GKbK 


Hk 
Gk 



{Bk - GK)bK \ 

-AKbK 

AkK 



BkK + ^ 





-2GKbK 
-Gk 
Hk 



(Bi 



In eqs. (B^-B^) we imply that x is the distance from the origin measured in 
units of m^; the energy eigenvalues E and the fermion mass mp are measured 
in units of m]y. 



For K = eqs. ( B.3 -B.9) reduce to those derived in ref. |]l6|; for any K but 
A = B = we get to the equations derived in a different context in ref. |31|. 

In order to find the Dirac eigenvalues numerically one can use various meth- 
ods. One method [^] is to express fermionic observables (such as the energy 
of the Dirac sea, etc.) through the phase shifts of the Hamiltonian TCk- In 
this work we shall use an alternative method: the diagonalization of the Hamil- 
tonian TIk ( |B.5| ) in the so-called Kahana-Ripka basis [^]. This method was 
previously used in the chiral nucleon model [^3|. The idea is to discretize the 
spectrum by putting the fermions into a large spherical box, and diagonalizing 
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the Hamiltonian in the basis of spherical Bessel functions which are eigenfunc- 
tions of the zero-field Hamiltonian. We will now introduce this technique. 

In the absence of the background field we have A(x) = 1, B{x) = 0, C{x) = 
0, H{x) = 1 and G{x) = 0, therefore Vk = 0, Wk = ^F, so that the free 
Hamiltonian becomes 

<) = ( ] , (B.IO) 

\^ mp -Vk j 

and the eigenvalue equation can be easily solved analytically. For given grand 
spin K and momentum p we find, generally speaking, eight linearly independent 
solutions of eq. ( |B.3p , which we call q = 1 ... 8; each of them is an 8- vector 
in the sense of eq. ( |B.4D . Introducing a column composed of the spherical Bessel 
functions, 

Jk.i = (jK+iipx), jxipx), jK{px), 3K-i{px), 

jx+iipx), jxipx), jxipx), jK~i{px)y i = 1...8, (B.ll) 

the eight solutions of the free eigenvalue equation can be written as 

v^^\ = MC^^^ Jk,! (no summation in i !) , (B.12) 



where M is the normalization factor (see below). With e = J + m?p, the 



(r) 

factors C\ are given by 



Cf' = ^ (-=^,0,0,0,1, ?,0,0|, 



c<.U-L(-o,„,„,:,_|,„.„) 

cp.-L(_E,_,o.„.„._!Ii..o,o).cf.-Lg,_,„.„.„-.„.„ 
cf.J.(„,„,-„.„,„,,_|), c«.-L(o.o.-!^.o,o,„.i.f). 
cf . -L („,„,_^._,„.o.„._!:^) . cf . -L (o,o.^,-, 0.0.0.^) ^ 

(B.13) 

It is convenient to discretize the free Dirac spectrum by the following trick 



[ p^ , 33 1 which preserves the spherical symmetry of the problem and makes the 

(r) 

eight states v)^ orthogonal to each other and to any other states with different 
K and e. Namely, we introduce a large radius R (eventually to be taken to 
infinity) and fix the spectrum for given K as zeros of the spherical Bessel 
function: 

miPuR) = 0. (B.14) 
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If 0(m,n are zeros oi Jk^z) one has the fohowing ortho-normaUzation conditions: 

/ dtfjK{amt)iK{ant) 
Jo 

= / dtt'^jK±l{amt)jK±l{ant) = 5mn7;[jK±l{oin)f . (B.15) 

Jo ^ 
These relations provide the orthogonality of the eight degenerate states 
( B.12| ) as well as their orthogonality to states with different p„ at given K. 
States with different values of K are orthogonal owing to the angular inte- 
gration. From eq. ( |B.15| ) we also learn that the normalization factor J\f in 
eq. pl2D is 

M = \jK±l{PnRT' . (B.16) 

We have thus constructed a complete ortho-normalized set of states uj^^ 
(given by eqs. dB.ll - BTl^ )) which are the eigenfunctions of the free Dirac op- 
erator. 

The eigenstates of the full Dirac Hamiltonian ( |B.5| ) can be now found from 
a direct diagonalization of the matrix 

'HK{Pm,Pn) = dxx'^vPi{pm,x)nK,ijVK]j{Pn,x) (B.17) 

(summation over i,j = 1...8 is assumed here). The superscripts r,s run over 



1 . . . 8 while the radial momenta Pm,n are given by eq. (B.14). They vary from 
small values of the order of 1/R up to some numerical cut-off Pmaxi which 
ensures the finiteness of the basis. Pmax and R have to be taken large enough 
so that no result changes when they are further increased. 



Appendix C 

In this appendix we derive matrix elements of scalar operators entering the 
Hamiltonians for bosonic and fermionic fluctuations about the spherically sym- 
metric sphaleron. These operators commute with the grand spin 

K = J-Ft = L-FS + t. (C.l) 

The fully normalized set of states in this momentum-adding scheme is 

\L,S,[J],T,K,K,) = Yl CfJ^!rnCiilssJSSs)\Tn)\LL,) . (C.2) 
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These states are the eigenfunctions of the operators K^, K3, L^, S^, and 
T^. On general grounds one can show that such Hamiltonians depend on the 
fohowing operators which commute with the grand spin: 

J2, V, t^, §2, S • L, S • T, t • n, S • n, T • [S x n], T • [n x L] , (C.3) 

where n = x/|x|. We use the general relations 

d % 

Any scalar operator commuting with the grand spin can be expressed through 



the set of operators ( C.3 ), for example 



t-L = ^(K2-j2-t2)-S-t or S-L = ^(J2-S2-L2). (C.5) 



To calculate matrix elements of the operators (|C.3| ) in the basis ( p.2| ) we use 
the technique of irreducible tensor operators and the Wigner-Eckart theorem 
(for details see Q). 

We illustrate this technique by calculating the complicated matrix element 
(L', 5', [J'], T', K, K3I T • (S X n) |L, S, [J], T, K,K^). The operator t does not 
act on the spin and angular variables so that the matrix element can be factor- 
ized as 

[J'],T', K, K3I t-(S xn) |L,5, [J],T,K,K^) 

= (-1)^+^'+^ I w' T } ^^'^'^'^'^ " ^ " " ^^^'^^^ " * " ' 

(C.6) 

where 



{T' II T II T) = ^T{T + 1){2T + 1) 5tt' ■ (C.7) 

The spin operator S acts only on spin variables whereas fi acts on angular 
variables, so one can again factorize the matrix element as 

{L'S'[J'] II S X n II LS[J]) = -iV2{L'S'[J'] \\ [S^^^ n(i)](i) || LS[J]) 

L' L 1 



i^6{2J + l){2J' + S' S 1 ^ (L' II n II L) (S' II S II 5). (C.^ 

J' J 1 
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Using the relation for the 9j symbol |34] 

L' L 1 
S S 1 
J' J 1 



{J' - L'){J' + L' + l)-{J-L){J + L + l) 
v/245(5 + l)(2S + l) 



one eventually finds 

{L', S', [J'],T', K, A'al t • (S X fi) \L, S, [J],T, K, K^) 

= i_^_^~^2J+2L'+T+S+l+K [(J/ _ j^i^^j, + ^/ + 1) _ (J _ + ^ + 1)] 



It{T + l)(2r + 1)(2J + 1)(2J' + 1)(2L + 1)(2L' + 1) 

J' J I \ { K t J' 1 ( L' I L 
LL'5|llJT|\0 00 



(C.IO) 



Along the same lines one can express matrix elements of operators of (C.3) in 
the basis (C.2) through 6j symbols; the result is: 



(L' , 5', \J'] ,T',K,K3\S-h\L,S,[J],T,K,K3) 



^_l)J+L+L'+s^jj^^(2L + 1)(2L' + l)S{S + 1){2S + 1) 

(C.ll) 





s 


v] 




• ') 




L 


« J 




y 



{L', S', [J'],T', K,K3\f-h\L,S, [J],T, K, K^) 

= |-_]^)2,/+2L'+5+T+ii'+l 



'{2L + 1)(2L' + 1)(2J + 1)(2J' + l)T{T + l)(2r + 1) 

\ K T J' \ { L' J' S \ I L' I L ^ ^ 





T 


J'] 




J' 


s' 








J 


TJ 




L 






y 



(L', 5', [/], T', K, K3I S • T |L, 5, [J], r, K3) 

= (^_iy+J'+L+S+T+K+l 



'{2J + 1)(2J' + 1)(2S + 1)S{S + 1)T{T + l)(2r + 1) 
■Sss'Stt'Sll' ,^r"i.^.r' (C.13) 





T 


^'1 


II 


Is 


J' 






J 


r J 




[■' 


5 
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{L',S',[J'],T',K,Ks\f-{hxl) \L,S, [J],T,K,Ks) 

= ]_^_iy+L'+L+T+S+K [^/(^/ ^ _ ^ -L)] 

2i 

■^T{T + l)(2r + 1)(2J + 1)(2J' + 1)(2L + 1)(2L' + 1) 

J' J 1 1 J K T J' 1 ( L' I L 
LL'5jjlJTj\^0 00 

We checked these expressions using the commutation relations: 

S-t,S-nj = -it-(Sxn), 

S-t,t-nj = iT-(Sxn), 

S • fi, t • nl = . 



(C.14) 

(C.15) 
(C.16) 
(C.17) 



For practical calculations the 6j symbols and Clebsh-Gordan coefficients can 
be evaluated using e.g. Mathematica. 



Appendix D 



In this appendix we state the divergent parts of the sea energy (4.1) in the basis 
of the set of eigenfunctions of the free Dirac-Hamiltonian, given by eqs. ( [B.llj - 
B.13| ). To this end a semiclassical expansion up to the quadratically divergent 
terms is performed. Subtraction of the result from the sea energy removes 
its quadratically divergent part for each value of the grand spin K separately. 
Since the quadratically divergent term is already complicated and lengthy in 
this basis, we do not continue the expansion to include also the logarithmic 
divergencies. Instead we use a simpler distribution of the total logarithmic 
divergence among the K sectors which may leave the separate sectors logarith- 
mically divergent, but renders the sum over all K finite. Since the total value 
of the logarithmic divergence is small, the complete removal of all divergencies 
for each K is of little use and would cause an enormous increase of numerical 
effort. 

We start from eq. (4^) , insert eqs. ( [B.12 , B.13 ), perform the sum J2r 
obtain: 



1 



div 



dt X - 

+3/2 



2e-*Pn 



R 8 



dxX^^JK,i{PnX) 



i=l 



JK,i{PnX] 



(D.l) 



div 
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We use the expansion 



(D.2) 



The constant 1 is cancelled by subtraction of the vacuum (free field). Con- 
tributions to the quadratic divergence arise from the terms 0{t), 0{t^p'^) and 
0{t^K'^). The result is as follows: 



with 



F{x) 



+0i - JoKA' - BC) + (-3i^ - jo^) - 



+(1 - 5k,o) 



2m^(ix+i + + jK-i){H' + 1). 
K 



f K + l .2 



\2K + 1 



K 



2K + 1- 



{A' - BC) 



{K + l){2K + 3) . . K{2K-1) , \ A' 



+ 2 ^ '^^.0 



+(1 - fe,o) 



{K + 1){2K'^ + 7K + 9) .2 



2K + 1 



j^+i + (2K^ + 2K + l)jl 



34 



+ I TTTT^— Jk+i Jk 



4{2K + 1) 
8K^ - 20K^ - UK + 5 .2 



C 



(D.3) 



A{2K + 1) 

where we have omitted the argument (pnx) of the Bessel-functions and denoted: 

X X X X ax 

To see whether this expression corresponds to the quadraticaUy divergent 
part of eq. (O) we perform the hmit i? ^ oo in which jK^i{pn)R^ R/Pn 



and l/-RX]n ^1/'''"/ '^P- Using 

Y,{2K + l)ii(px) = 1, Y.^2K + Ifjlipx) = 1 + \{pxf, 

K K 
and performing the p integration we obtain: 



(D.5) 





dt 




vr Jo 




/ 






dt 








47r J 


'0 


T 



dxiA^ + B^ + 



5C2 



(D.6) 



The first term coincides with the quadraticaUy divergent part of eq. ([4.3D , as it 
should do. The second term is logarithmicahy divergent so that the remaining 
logarithmic divergence of eq. ( |4.3| ) to be removed is therefore: 

Ediv,i = 77- / "T / dxg{x), 
Zn Jo t Jo 



gix) 



X \ X Ix'^ 



+'mpx'^ 



{G' + H'^f - 1 



mi 



+2x2 (G'^ + H'") - 2xC{HG' - H'G) 



(1-^)2 + B2 + 5c! 



2x2 



(D.7) 



This divergence can be distributed among the K sectors as follows: 

1 ^ p-*Pn rR 



Edi.,i = ^(2i^ + l) . 



dtt^'^Y. 



\jlk+l{PnR)B? Jo 



dxx G{x) 
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with 



G{x) = 6k,0 [jliPnX] +jo{PnX)^ 



+ (1 - Sk,o) {jK+liPnX) + 2i^(p„x) + ^ . (Di 



As already mentioned, we would obtain a much more complicated distribution 
of the logarithmic divergence among the K sectors if we continued the expansion 
eq. ( D.2 ) up to fourth order and collected the logarithmically divergent terms. 
Although the simpler version eq. ( p.SD removes the logarithmic divergencies 
only after summing over K and not for each K separately, it is sufficient for all 
numerical purposes to use this formula. 

Let us finally remark that one can in principle use a simpler formula also for 



the quadratically divergent terms which can be obtained by expanding eq. (p. 2 



in a large K limit. This leads to an expression for Ediv,2 similar to eq. ( p.3[ ) 
but contains only the first term + — 1. Numerically identical results 
are obtained with this version, but one needs higher numerical parameters to 
insure stability, so that the numerical effort is increased. 
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Figure captions 

1. Energy barrier as a function of Ncs for different values of the chemical 
potential fj, for niH = my/- The units for EdassifJ-) and /i are Mq = ficrit = 
2TTm\y /a. 

2. The discrete level e^ai and some discretized continuum eigenstates as a 
function of Ncs for rnp = mw- It is demonstrated that each level is 
shifted upwards and finally takes the position of its predecessor. 

3. The renormalized fermionic energy Ej"^ as a function of Ncs for the 
fermion masses mF/miy = 0, 1 and 2. It has to be compared with the 
classical energy which is about 100 niw for Ncs = 0.5 . 

4. The fermionic temperature dependent part PE^^^ (solid lines) and the 
classical part PE^f^^^ (dashed lines) of the sphaleron transition rate as a 
function of T for rrit/mw = 1-5, 2.0 and 2.5. E^f^^^ and E^^^^ are defined 
in eqs. (g, |]2l|). 

5. The baryon density as a function of the radial distance r for mp = rnw 
and Ncs = 0.26, 0.5, 0.74, and close to 1. 
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Minimal Energy Paths for Different /i 
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